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ABSTRACT 


The image reconstruc tion techniques used widely in 
the areas of medical imaging has been successfully 
adapted to measure the void fraction and attenuation 
constant of the material. The convolution backpro jection 
algorithm in conduction with the Ramchandran~ 
Laxminarayanan filter has been used for reconstruction 
in the form of a density distribution of a 
cross-section. The observation of the reconstructed 
results shows that the process of reconstruction is more 
vulnerable to uncertainties inthe region close to the 
centre of the objects. The uncertainties arises due to 
the error incorporated in the projection data and the 
weaker sources. 


4 



ACKNOWLEDSEMENT 


I wish to express my indebtedness to Dr.Prabhat Munshi for his 
invaluable guidance, his accessibility and unlimited co-operation 
during the course of study. It was due to his moral support that an 
uninitiative like me could attempt a project in this fascinating 
field.lt was a nice experience to work under him. 

I sincerely acknowledge the support of the entire NET family, 
especially Mr . S . S .Pathak for his invaluable guidance. 

I am grateful to my family for consistently encouraging me 
despite the distance. 

I can not forget the tasty meals I had, thanks to Mrs. Munshi. It 
certainly alleviated the monotony of the hostel food. 

My sincere thanks to my friends and colleagues for the helpful 
discussions, inspiration and co-operation throughout ray stay over 
here . 


Q . K . Dwivedi . 



TABLE OF CONTENTS 


Abstract 4 

Acknowledgement 5 

List of figures 8 

List of symbols 9 

1.0. Introduction 10 

2.0. Preliminaries 12 

2.1. Data collection modes 13 

2.2. Parallel beam geometry 14 

2.3. Fan beam geometry 16 

2.4. Image inversion formula 18 

2.4.1. Parallel beam formulation 18 

2.4.2. Fan beam formulation 20 

2.5. Convolution backpro jection algorithm 22 

2.5.1. Parallel beam algorithm 22 

2.6. Radiation statistics 24 

3.0. Experiment using single channel analyzer 27 

3.1. Equipments needed 27 

3.2. Description of equipments 27 

3.2.1. Scintillation detector 27 

3.2.2. Linear amplifier 31 

3.2.3. High voltage unit 31 

3.2.4. Low voltage unit 31 

3.2.5. Single channel analyzer 31 

3.2.6. Timer 32 

3.2.7. Scaler 32 

3.2.8. Gamma-ray source 32 


6 



4.0. Program implementation and description 34 

4.1. Computer implementation of the CBP algorithm 34 

4.2. Program for reconstruction of the image from data 36 

5.0. Results 39 

5.1. Data used 39 

5.2. Results 40 

6.0. Conclusion and recommendation 47 

References 51 

Appendix 

A. Data employed 53 

B. Pictures of objects 57 

C. Computer program for image reconstruction 62 


7 



LIST OF FIGURES 


1 . Parallel beam geometry 
2 • Fan beam geometry 

3 . Conversion of fan beam geometry 
to parallel beam geometry 

4 . Reconstructed image of object 1 , scan 1 

5 . Reconstructed image of object 1 , scan 2 

6 . Reconstructed image of object 2 , scan 1 

7 . Reconstructed image of object 2 , scan 2 

8 . Reconstructed image of object 3 , scan 1 

9 . Reconstructed image of object 3 , scan 2 


15 

17 

23 

41 

42 

43 

44 

45 

46 


8 



list of symbols 


A 

B 


f(rA) 

f(r»#) 

1 

M 

N 

No 

p( s,®) 
r ,♦ 
s 
s' 

As 

W(F) 

0 


Fourier frequency; 

Angle of extreme data-rays; 

Chord of integration; 

Distance of source from origin; 
Partial derivative of '» 

Function being reconstructed; 
Reconstruction of /(r,#)i 

projection data for fan beams; 
Variable of integration; 

Thickness of absorbing medium; 
Radiation detector reading; 
Radiation source strength; 
Projection data for parallel beam 


Cylindrical co-ordinates; 
perpendicular distance of data 
s-distance of data-ray passing 


from origin; 
through (r,4) 



spacing between 


data-rays; 


Window functxoH) 

Source position(FAG ) y 
source position(PEG) ; 
absorption co-efficient. 


9 



CHAPTER 1 
INTRODUCTION 


Image processing has revolutionized the area of medical 
imaging, which is nowadays used in a variety of ways in non 
destructive testing(NDT) , multi phase flow. Nuclear engineers have 
used image processing techniques in the measurement of void 
fraction during "LOSS OF COOLANT ACCIDENT". Because they are very 
reliable & accurate. 

The experiments, which are performed by us using single 
channel analyzer ( SCA) , we measured cross-sectional distribution of 
the objects indirectly by probing the objects with the collimated 
invisible radiations and then interpret these results. We have 
related the measured data to the cross- sectional distribution in 
a known way. In the image reconstruction procedure, we processed 
the data to form the cross-sectional image i.e. interpreted the 
results indirectly. We take various strip integrals corresponding 
to a particular angle of view. This set of integrals is called a 
projection of the object. Given a number of such projections at 
different angles of view, the estimation of the corresponding 
distribution within the object is the basic problem of image 
reconstruction from projections. 
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The principles, which are used to reconstruct the image, are 
attributed to Radon[2],who proved that any arbitrary function could 
be recovered from its set of line integrals taken along 
various chords & directions. However, the method of Radon could not 
be implemented due to inherent mathematical complexities . Bracewell 
reported an application in radio astronomy and CormackCS] derived an 
inversion formula, both of which were closer to being implemented as 
compared to Radon’s solution. Bracewell & Riddle! 12] and Ramchandran 
& Laxminar a yan[ 9 ] showed how computations involved could be 
accelerated by adopting the use of convolutions. 

The image processing can be used in diverse & wide ranging 
areas like astronomy & electron microscopy. We can use image 
scanners to estimate the quality of meat prior to slaughter. It 
can also be used to inspect wooden poles in power transmission. 

The present study is an effort towards studying the image 
reconstruction through convolution backpro jectlon algorithm which 
can be extended to various areas of applications. 
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CHAPTER 2 
PRELIMINARIEB 


Single beam mono-energetic radiation, attenuation phenomenon in 
a plane can be represented by 

N = Nq exp[-f^Ji(r,^)dll (1) 

The value of H is characteristic of a material for a given type of 
radiation and the energy of radiation. Equation (1) considers M to 
be a two-dimensional function of positioni as the path of radiation 
is assumed to be restricted to a plane. The energy dependence of 
is also ignored. The simplifying assumption is valid for 
mono energetic radiation sources, e.g. y-rays, and suitable 
corrections are incorporated to handle X-ray sources which have a 
distributed energy spectrum. We note that when the cross- section 
of interest has a uniform distribution of the absorbing material 
then (1) reduces to, 

N = Nq exp ( -mAI ) , ( 2 ) 

where g is no longer a function of position. 
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In general, any cross-section of interest will have non- 


uniform distribution. Rewriting (1) as, 


p = J /i(r,#)dl 


(3) 


where 


p = ln(NQ/N) 


We observe that recovering M from the given data p is not a 

straightforward provision. Radon showed that it is possible to recover M 
from a set of several p-values measured along various chords, c. 

These /n-values can then be suitably calibrated to give the density 

values if so desired. 


2.1. DATA COLLECTION MODES: The image processing methodology requires 
attenuation data collected by an array of radiation detectors for 
the reconstruction of the function Two popular data 

collection modes are termed as parallel-beam geometry( PEG ) and fan 
beam geometry ( FBG) . 
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2. 2. PARALLEL BEAM GEOMETRY; The IP scanners based on PBG mode 


have several pairs of radiation source and radiation detector 
systems which scan the object completely. Figure 1 depicts the 
configuration for a generic scanner. The source-detector pairs are 
spaced uniformly and the object to be imaged is stationed on a 
table which can be rotated to give different values of The 

source-detector pairs are spaced uniformly and the object to be 
imagined is stationed on a table which can be rotated to give 
different values of The line SD represents the path of the data 
ray. The perpendicular distance of the ray from the center of the 
object, which happens to be the origin of the frame of reference 
chosen for derivational and computational convenience, is denoted 
by s. Several SD pairs collect the data p for a given This set 

of p is known as 'projection*. The object table is rotated to get 

several sets of p for different values of 0, i.e., different 

’views*. The data thus collected is denoted by p(S|;0),the 
semi-colon indicating the procedure of data collection. If the 
object is stationary, then the SD system has to be rotated to 
collect the data for various ’views’. 
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S - Source 
D -Detector 


RAY = (s , 8) 
f = f (r,0 


F'IQ. 1 Porollel beam collection geometry. 


2. 3, FAN BEAM GEOMETRY ; The medical scanners incorporate the FBQ 
collection mechanism which involves a rotating module including a 
source and an array of detectors stationed in an arc around the 
patient. Figure 2 shows the FBG configuration. Here, a single 
source views several detectors simultaneously and data are 
collected by different detectors corresponding to various angles, 
flf.This process is repeated for different values of P to to get the 
IP data, g(ff,/J).The FBG mode can also be implemented for a 
rotating object system also. The data collection mechanism remains 
the same as that for the PBQ case. 

2. 4. IMAGE INVERSION FORMULA; The process of reconstructing a 
two-dimensional function from its projection data is similar for 
the FBG and PBG cases. The parallel beam IP formulation involves 
Fourier transforms while the fan-beam formulation incorporates 
Hilbert transforms in the inversion process. 


2.4.1. PARALLEL-BEAM FORMULATION: The inversion formula for the. PBG 
case is based on the ’central-slice* theorem for Fourier transforms 
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Fig. 2 Fon-beam data collection geometry. 
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of two-dimensional functions. The theorem states that the 
one-dimensional Fourier transform of the projection data, 

p(s;®), with respect to the first variable s, is equal to the 
two-dimensional Fourier transform of the object function f being 
imaged. Mathematically, 

f(F;^) = p(F;0) . (5) 

talcing the inverse fourier transform of (5) we get 

f(r,#) = fo^-oo P( F»®)exp[-i2nFr cos (0-^)]|FtdFd6. (6) 

Here for a given 0, 


P(F;^) = f^j^p(s;6)exp[-i2itFs]ds (7) 

Equation (6) requires projection data on a continuous basis for 
all values of s and ®. For computational feasibility a finite 
cut-off is implemented in the form of a filter function, 
W(P), which vanishes for |fI greater than A, the cut-off 
frequency . Thus the filtered version of (6) becomes, 
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{ 8 ) 


f{r,#) = JqJ *«5 p( F;®)exp[”i2npr coa(0-^)]W( F)jFjdFd0. 

Here, the reconstruction, f, is approximate due to the finite 
cut-off introduced in the computations. Also, the discrete nature 
of the IP data forces the practical Implementation of A to be 
governed by the sampling theorem, i.e., 


A> [1/2As] 


The band-limiting filter, introduced in electron micrography by 
Ramachandran & Laxminarayan is given by 


W(F) s ( 


l,iFl<A, 

0,lF|>A, 


and a popular ’sine’ filter used by Shepp & Logan is given by. 


W(F)k 


^[sin(IlF/2A)]/(nF/2A),|Fl<A 

(0 , 1f1>a 


2. 4, 2. FAN BEAM FORMULATION: The Inversion formula for the FBG case 
was first derived by Herman & Naparstek and is given by, 
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/(r,*) = (l/4ir2)f2njB^tl/sin(^..^)D^^(or,p)d(Td^, 
where , 

j?(®'>/*)=data for the ray represented by 
T)^g{<^J) = {l/U)[{9gAr)-{9g/d^)], 

(»■’= tan [rcos(P-#) ]/[D+rsinO“^) ] } , 

U ={ [rcos(^-^) ] [D+rsin(^~#) , 

and other variables are as in the PBG case. We note that U is the 
distanse of the radiation source from (r,^),the point being 
reconstructed* and la the angular displacement of the 

particular data ray passing through that point (r*#).A major 
difference is the computation of partial derivatives of the 
data, g. 
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2 , 5. CON VOLUTION BACK PROJECTION ALGORITHM: The reconstruction 
algorithm based on the inversion formulae i are known as ’transform 
methods’. It is observed that these transform methods are slow and 
required very accurate interpolation schemes in computing the 
two-dimensional inverse fourier transforms. The introduction of 
convolutions eliminated the steps of computing the fourier 
transform of the data and the subsequent time consuming 
two-dimensional fourier inversion. 

2 . 5 . 1 . PARALLEL BEAM ALGORITHM: Exploiting the convolution property 
of fourier transforms ,{ 8 ) becomes as, 

^(r,^)*jQj\p(s;0)q(s’_s)d3de (10) 

where , 


q(s)=J^^W(F)lFlexp[i2nFs]dF, 


( 11 ) 


and, 


B* =rcos («_♦) 



Here, q known as the convolving function and is the inverse 
fourier transform of the filter function, W(F). Equations( 10 ) -( 11 ) 
are the fundamental equations of the convolution 

back-pro jection( GBP) algorithm of IP. The inner integral of (10) 
is a convolution and the outer integral was termed back projection 
by the early researchers. The GBP algorithm is the moat widely used 
method of reconstruction. 

2 . 6 .RADIATION STATISTICS :The quality of reconstruction is affected 
with four problems; 

l:Due to errors in data collection process, 

2:Due to the finite number of projections i.e. , discretization of 
the problem, 

3: Due to errors by numerical computations, 

4: Due to statistical nature of photons, because of their interaction 
with the matter and the photon detector. 
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In gamma-ray image processing, the rays are generally 
monochromatic or at most, two wavelengths are presents. Due to 
this nature of gamma-ray, uneven detector response to the 
detector, and of beam hardening, which affects the X-ray scanners, 
are resolved. 

If we consider that all the photons emitted by the source in a 
unit period of time in the direction of the detector are detected, 
the possible counts give rise to a discrete random variable 
denoted by M. It can be further shown that the probability of 
occurrence of M.It can be further shown that the probability of 
occurrence of M, at a specific value m is denoted by, 

P^(m)=exp(-J:)x"*/m! (12) 

where ^ is a fixed real number. 

Equation (12) is called the poisson probability law, and M which 
satisfies this law is called a poisson random variable with 
parameter The main properties of this variable are as follows: 

(i) its mean is X, 

(ii) its standard deviation is Vjf, 

(iii) it behaves normally if 71. is (>100), 
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X. will be estimated by the count of the the number of photons 
during a particular period of unit time i.e., by a sample of the 
random variable. Poisson statistics imply that if we increase the 
sample by a factor of N, we increase the sample by a factor of N 
we reduce the size of the l<f error by a factor of , where is 
the standard deviation. 

A photon leaving the source in the direction of the detector, 
will reach the detector (without being absorbed or scattered) 
with a fixed probability .This probability depends on the energy 
of the photons and the material lying on the line between the 
source and the detector. is called the transmittance of the 
material (along that particular line), at that particular energy 
i.e., the photons which leave the source in the direction of the 
detector, a fraction will eventually reach the detector, the 
rest being absorbed or scattered. The photons on reaching the 
detector are counted with an efficiency t, which is called the 
efficiency of the detector. 

Thus the number of photons which actually reach the detector 
without being absorbed or scattered, and are counted by the 
poisson variable with parameter 


25 



As discussed earlier, 


„ln(N/NQ)=f^M(r,^)dl s 

where S is the ray sum & is the sample of a random variable such 
that , 


|Mg+ln(N/NQ)|< T, 

Where T =* <Va^a*d^“^ 

and is the fraction of photons leaving in the direction of the 
detector , 

% is the number of photons emitted during the period of 
d 

measurement , 

^ is the transmittance of the material and 

Ctt 

is the efficiency of the detector, 
d 
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CMAPTER 3 

EXPERDBIT USING SINGIE CHANNEL ANALIiZER 


3.1. The equipments needed tor the experimental set up are 

1. Nal scintillation detector 

2. Linear amplifier 

3. High voltage unit 

4. Low voltage unit 

5. Single channel analyzer 
6 .Timer 
7 . Scaler 

8. Gamma- ray source 


3 . 2 .Description of equipments; 


3. 2 . 1 . SCINTILLATION detector: There are certain materials 

(e.g.,organic.inorganio,and plastio)«hich emit flashes of light or 

' ♦••inn -oasses through them. These 

aolntlllatlons when radiation p 

scintillators can be used to detect any type of radiation hut are 

most Widely used for gamma-ray counting. The output of a 

4 . • «oi +n the incident energy>and hence can 
scintillator is proportional to the incia 

be used to identify the energy of a gamma-ray. 
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Radiation is absorbed in the scintillator when it strikes the 

material . The absorbed energy is utilized in lifting the electrons 

to the excited state and part of it appears as heat energy. The 

electrons in the excited state de-excite in a very short timei 

leading to the emission of light. The de-excitation time or the 

time constant of a scintillator varies with the material land the 

—8 

decay time for commonly used scintillators is of the order of 10 , 

The intensity of the emitted light is directly proportional to 
the energy lost by the incident particle. The emitted light flashes 
are amplified by a photo multiplier . The photo multiplier tube 
consist of a photo cathode and several dynodes, each of which is 
maintained at a successively higher voltage. The photo cathode is 
usually a thin coating of antimony-cesium on the inside surface of 
the photo multiplier , whereas the dynodes are 

antimony-cesium, silver-magnesium, or copper-beryllium, The electrons 
emitted by the photo cathode are amplified by releasing more 
electrons at each of the dynodes, The electrons from the last dynode 
are focused on to the the anode connected to the preamplifier. 
The coupling of the scintillator to the photo multiplier is done by 
a light pipe. Normally, the scintillator is covered with a highly 
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reflecting material, such as MgO and A1203, to avoid loss of light 
pulses. The light pipe is commonly made of glass, Incite , or 
plexiglas. Its main purpose is to transmit efficiently the light 
pulses to the photo cathode of the photo multiplier tube. 

The advantages of using a scintillator for detecting 
gamma-rays are as follows: 

1: Since all the light pulses are emitted in a very short time, 
the scintillator integrates these into a single pulse, The height 
of this pulse is proportional to the incident gamma-ray energy. 

2: Since the scintillator is solid, it has more atoms per unit 
volume than a gas-filled detector, leading to a better detection. 

3: Since the time constant . of the scintillator is small, its 
dead-time is very small, facilitating a high count rate. 

4: A large scintillator increases the probability of absorbing 
the scattered gamma-rays, leading to an increase in photo peak 
efficiency. 
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A disadvantage in using a scintillator is the decrease in the 
light output with increase in the energy of the incident 
gamma-radiation. For example, the Nal (T1 ) scintillator displays 20 
percent variation in the light output per unit energy of gamma-ray 
over the range 0-1 MeV , The light output also varies greatly with 
the type of radiation. With the Nal(Tl) scintillator, electrons 
give a better response than alpha-particles for the same energy. 

The overall efficiency of a scintillator is low because the 
detection process, namely, radiation cause light emission which 
is then reconverted into electrons at the photo cathode, is 
complex. Typically 300eV of energy is dissipated in the NaI{Tl) 
scintillator for the release of one photo electron; in an organic 
scintillator, the energy required is l-5keV per photo electron. The 
corresponding energy for producing an ion-pair in a gas- filled 
detector is 30eV and in a semiconductor is 4eV. 

Scintillators are many types, e.g., disk-type, right circular 
cylinder, and well-type. 
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3. 2. 2. LINEAR AMPLIFIER: It is an amplifier, which has three 
gain controls. 

(i) Input attenuator, which have attenuation factors 1, 2, 5, 10 

20, &50. 

(ii) Coarse gain, which have amplification factors 0.1, 0.3, 1. 

(iii) Fine gain, with an amplification 1-3. The pulse shaping 
provision of linear amplifier has variable integrating time 
constants O.lM sec, 0.2fi sec, 0.5ft sec, Ift sec, 2/t sec& 5|t 
sec. Also it gives output pulse of amplitude 8V. 

3. 2. 3. HIGH VOLTAGE UNIT: It is a high voltage supply of 300 to 2500V 

3. 2. 4. LOW VOLTAGE UNIT: It is a low voltage supply of _24 to +24V. 

3. 2. 6. SINGLE CHANNEL ANALYZER ( SCA) : An SCA can operate either in the 
integrating ( INT) or differianting(DIFF) mode. In the INT. mode, 
the lower level (base line) of the SCA acts as the discriminating 
bias. In the DIFF mode, the lower level is set at E & the upper 
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level at E+Ae so that only the signal with a pulse height between 
E & E+AEis counted by. the scaler. The window Ae can be 
continuously varied from 0 to 2V. The base line adjustable between 
0.2 to lOVjThus making it possible to resolve & identify any pulse 
amplitude in this range, the output of SGA is constant & of 
magnitude lOv. 

3 . 2 . 6 .TIMER: It is an ordinary counting set up to count the time 
required to take the counts. The factors are 10 sec> 50 sec, 100 
sec, 500 sec, lOOOsec. 

3.2.7. SCALER: It is a counting set up to count the signals of pulse 
height between E & E+Ae. 


. , . 60 ^ 

3 . 2 .8 .GAMMA- RAY SOURCE: The y-ray source used by us is Co. 

( iii) PROCEDURE TO COLLECT THE DATA:To collect the data we have to 
go through two steps: 

1 : IDENTIFICATION OF THE PHOTO PEAK OF A RADIOACTIVE SOURCE: Connect 
the low voltage and high voltage unit to the preamplifier of the 
scintillation head. Connect the output signal of the preamplifier 


32 



to the linear amplifier .Connect the amplifier output to the single 

channel analyzer { SCA ) . Connect the SCA output to the scalar 

through the timer. Turn on the low voltage unit, wait for a few 

minutes, and turn on the high voltage unit. Slowly increase the 

detector voltage to the value 500V. Place the source through a 

collimator. Set the window Ae to 0 . IV & slowly increase the base 

line from O.lV.At each base line setting, for 10 sec time, collect 
the data The counts will be high near the photo peak of the source. 

2. COLLECTION OF THE DATA AT THE IDENTIFIED PHOTO PEAK; Set SCA at 
the voltage i.e.the voltage at which peak is identified of V~ray 
source ^Co. Using the experimental set up as in the stepl , 

collect the data using the object under parallel beam data 
collection mode. 
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CHAPTER 4 


PROGRAM IMPLEMENTATION AND DESOHPTION 


There are various programs developed, which simulates the 
projected data, reconstruct images from actual or from simulated 
data. These programs have been written in Fortran. 

4.1. COMPUTER IMPLEMENTATION OP THE CBP ALGORITHM: 

The problem of reconstruction from projections is as follows: 

Given p(s;<?) and, find f’(x,y) i.e., given discrete 

projection data in the form of the estimates of p for a finite 
number of rays; find a 2-d distribution , which is a reconstructed 
estimates of the unknown object. 

In the case where p is sampled uniformly in both s and 0, for 
N angles apart, with each view having M equispaced rays As 
apart, we define 

M‘‘‘= (M-1)/2 

}M odd 

M"=-(M-1)/2 
M‘^= (M/2)-1 


M =-M/2 


}M even 



In order to ensure that the collection of rays specified by, 


{(mAs, nAe):M"s <n>N } 

covers the unit circles, we have, 

A0=*/n and As = l/M^ 

A reconstruction algorithm which can be implemented on the digital 
computer is required to evaluate f(kAx,lAy), which is a band 
limited approximation of the function to be reconstructed. Here 
K”s k 2:K^& 1 & L~ , where k and 1 are the position of the 

co-ordinates of the image pixel. The definition of their upper and 
lower limits is similar to that of . Thus the projected data 
from N views 7 M rays is to be used to construct an image of KxL 
pixels. In this particular case, the image is composed of 21x21 
pixels from 18 views each having 21 rays. The back projection 
integral is evaluated as follows: 

N 

f(kAx,lAy)2‘ 3^ p(kAxoosfi +1 Ay sin®) 
n=l 

For each angle ® the convolved values of p(s,® ) for the 

n , 

KxL Values of s. We can either have a separate convolution for 
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every s’ with the actual value of q(s’-mAs) at that point or we 
can evaluate p(tnAs, 0^) only within the specified limits of m & 
then use interpolation . The latter approach is much faster and 
cheaper. These operations are represented by, 

M ^ - + 

p (mAs,6 As X p(raAs,6n) q(m-inAs), M 5 m > M 

” + 
m=M 

Pj.{s’,0j^)^ As S p^(mAs, 0^) I(s-roAs) 
m’ 

where I(s) is an interpolating function, 
interpolating function, say I(s), corresponding 
interpolation between adjunct samples is 

I^vs)- ‘^1/As( 1-lsl/As ) Isis As 

using the above formulae ^ a program was written in FORTRAN to 
implement the CBP algorithm# 

4*2. PROGRAM FOR RECONSTRUCTING THE IMAGE FROM DATA: 

The program (see appendix B) can reconstruct an image from 
experimental data. The program is written in Fortran to accept the 


A linear 
to linear 
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da'ta. for the parallel beam geometry and use to reconstruct the 
image. Images of several experiments are reconstructed within 
acceptable error limits. 


A suitable image quality index would be helpful in 
quantitatively estimating the accuracy of image reconstruction. Two 
such indices are l^and Ig, The process of reconstruction can be 
terminated if the error of reconstruction reaches a certain minimum 
value. The actual program details are follows: 

1. The parameters of the program LITM(the number of rays per 
projection) and LITP( parallel beam data) are read in from data 
file , 

2. The projected data are read in from a data file. This 
transformation enables us to use the CBP algorithm* for parallel 
beam data* without resorting to the Hilbert transform. 

3. The Ramchandran- Lakshminarayan filter is used in this program. 
The numerical values had been computed and stored in a file. These 
values were merely read from the file into an array. It is at this 
stage that the convolution is carried out. 
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Convolution is especially easy if the values to be convolved 
are stored in an array. It is often found that, during 
back-projection of data, the the convolved values do not exactly 
pass through the point being reconstructed. Due to this 
interpolation has to be used. Since linear interpolation is fast 
and it gives good results, it has been used in this particular 
implementation • 

4. The superimposition of these interpolated values is carried 
out. Some convenience like display of point values as they are 
reconstructed have been added to the program. These reconstructed 
values are the LITF values. 

A sample output of this program has been displayed in appendix 

C. 
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CHAPTER 5 


RESULTS 

In this chapter, the results obtained with the actual 

experimental data have been presented. 

5.1. DATA USED: The data used are taken from the study conducted 
by us using single channel analyzer( SCA) . The details are 

briefly summarized as follows: 

a) The source of nuclear radiation is 1.33 Mev of Co-60. The 

collimator is fixed at the other end. The source and detector are 
placed at the fixed place but the object is moved, perpendicularly 
to the axis of the collimator, which covers the whole outer 

diameter of the object, because of parallel beam formulation. 

b) The detector is connected to the other end of the collimator 
and the subsequent stages being an amplifier, single channel 
analyzer, a timer and a scaler. 

c) The settings of the apparatus are: 

Polarity is = -ve, 
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Attenuator = 20, 

Gain = 0,16, 

High voltage = lOOOV, 

Lower level of window = 0.5V, 

Window width = O.IV. 

d) The data collected for calibration are for a cylinder, a small 

cylinder in a large cylinder, and water in beaker, Raw data appears 
in Appendix A. 

5.2. RESULTS: The reconstruction results are given in Fig 4-9. 

The results shows that we can reconstruct the image of the 
object of any shape and sort out the defects in the object through 
nondestructive testing. We calculate the void fraction during a 
"LOSS OF COOLANT ACCIDENT" ( LOCA) . We can measure 
cross-sectional distribution of any property through indirect 
measurement by pi'obing the object with invisible, penetrating 
radiation and interpretation of these results. 

Expression (13) in section 2.6 shows that the error in 
measurements depends on various different factors and it may be 
reduced by , 
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NEGATIVE IMAGE 

image datal (hallow cylinder) 
ra(ii257 
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Fig. 4 
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NEWTIVE IMAGE 

image dataH (hallaw cylinder) 
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Fig, 5 
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NEGATIVE IMAGE 

image data3 (water in a beaker) 
ra(i\S57 
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Fig. 6 
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image data 4 (water in beaker) 
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Fig. 7 
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NEGATIVE IMAGE 

image data 5 (small cylinder in a large cylinder) 
rani257 
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Fig. 8 
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NEGATIVE IMAGE 
image data 6(siiiali cylinder in a large cylinder) 
ram257 
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11 

11 

11 

11 

11 

11 

11 

11 

11 

11 

36 

36 

38 

32 

29 

20 

18 

19 

33 

11 

11 

11 

11 

11 

11 

11 

11 

11 

11 

11 

11 

11 

11 

11 

11 

43 

11 

11 

11 

11 

11 

11 

11 

11 

11 

11 


MINIMUM LITF=-0.61044890E-01HAXIKUH LITF= 0.48694250E+00 
HORIZONTAL CENTERLINE 

\ 

1 11 0.14145230E+00 

2 11 0.14603000E+00 

3 11 0.23134540E+00 

4 11 0.2829Q670E40Q 

5 11 0.26219300E+00 

6 11 0.25808520E+00 

7 11 0.36757A40E+00 

8 11 0.24262980E+00 

9 11 0.23847800E+00 

10 11 0.10640E00E+00 

11 11 0,486942502+00 

12 11 0.98800910E-01 

13 11 0.22883360E+00 

14 11 0.24277430E+00 

15 11 0.38989580E+00 

16 11 0.25710580E+00 

17 11 0.28732000E+00 

18 11 0.25873290E+00 

19 11 0,21273050E+00 , 

20 11 0.97292160E-01 

21 11 0.174996aOE+00 


AVERAGE = 0,228048702+00 



1. Increasing the strength of the source or the time period of 
measurements. 

2) Improving the efficiency of the detector. 

The projection for the objects are assumed to be less affected 
by statistical errors since averaging over a large number of 
readings has been carried out. 

The detector for photon detection has been collimated to 
accept only those photons which come straight from the source. 
However, perfect collimation is difficult to achieve. This allows 
some scattered radiation to enter the collimator window, 
artificially boosting the counts. Since only straight attenuation 
photons are related to the chordial density, scattered photons 
produce an overestimate of the density(LITF) . These effects due 

to photon scattering have not been taken into account. The 

increased number of counts due to scatter will not results in an 
overestimate of if the effect of scattered photons is the 

along all chordial directions. This is a fairly strong assumption 
and may not be always valid. It has possibly resulted in 
underestimates of density (LITF), partially compensating for the 
greater effects of counting statistics at lower counts. 
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CHAPTER 6 

CONaCSION AND RECOMMENDATION 


It has been observed in the course of study that image 
reconstructed can be extended to get the meaningful results for 
measurements of void fraction in a multiphase flow. It is also 
possible to have LITF to be identical to the attenuation 

constant of the material at that energy, by suitably normalizing 
the input data and ensuring that the same system of units is 
consistent during computation. 

The error in reconstruction occurs due to the error 
incorporated in the projection data, and the error in the 
projection data occurs due to counting of scattered 
photons . 

The results also shows that the error occurs due to imperfect 
collimation of the source and the poor activity of the 
source. The accessibility of the experiment can also be extended by 
use of automatic system, which can be done by use of feedback 
using PC. 

It is apparent from the above recommendations that further 
investigations are needed along these lines. Besides, the study 
is conducted using only Ramchandran™ Lakshrainarayan filter. 
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Investigation for other window functions can be carried out using 
limited amount of data. A larger number of rays will also give 
a better results. Of course this work has been done earlier using 
automatic system but whatever we have done using manual system 
gives good results as comparable to the automatic system. 
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APPENDIX A 

Lower level voltage 
0.0 
0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 

1.0 

1.1 

1.2 

1.3 

1.4 

1.5 

1.6 

1.7 

1.8 
1.9 
2.0 
2.1 

2.2 

2.3 

2.4 

2.5 

2.6 

2.7 


Counts per 10 secs 
76 
30 
27 
27 
26 

25 

26 
29 
29 
39 

47 

50 
53 
53 

51 
68 
61 
53 
51 
49 

48 
48 

51 
66 

72 (photo peak) 
58 

52 
47 
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OBJECT 1. 
HOLLOW CYLINDER 


DISTANCE MOVED BY OBJECT ( eras . ) 

AIR 

0.0 

0.5 

1.0 

1.5 

2.0 

2.5 

3.0 

3.5 

4.0 

4.6 

5.0 

5.5 

6 . 0 

6.5 

7.0 

7.5 

8.0 

8.5 
9.0 

9.5 

10.0 


COUNTS PER 50 SECS 


SCAN 1 
1637 
1545 
702 
874 
977 
1118 
1124 
1134 
1179 
1154 
1205 
1114 
1184 
1146 
1152 
1070 
986 
994 
788 
1065 
1486 
1602. 


SCAN 2. 
1648 
1532 
745 
893 
950 
1123 
1153 
1143 
1198 
1165 
1234 
1136 
1167 
1178 
1186 
1076 
974 
1004 
796 
1076 
1478 
1611 
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OBJECT 2 


WATER IN A BEAKER 

DISTANCE MOVED BY OBJECT(cms.) COUNTS PER 50 SECS. 



SCAN 1 

SCAN 2. 

AIR 

951 

941 

0.0 

814 

851 

0.4 

727 

687 

0.8 

656 

642 

1.2 

627 

572 

1.6 

628 

596 

2.0 

620 

591 

2.4 

573 

554 

2.8 

565 

557 

3.2 

548 

528 

3.6 

575 

556 

4.0 

562 

586 

4.4 

594 

566 

4.8 

568 

577 

5.2 

560 

627 

5.6 

539 

565 

6.0 

609 

666 

6.4 

620 

671 

6.8 

715 

775 

7.2 

808 

879 

7.6 

890 

84 9 

8.0 

812. 

852 
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OBJECT 3 


SMALL CYLINDER IN A HOLLOW CYLINDER 


DISTANCE MOVED BY THE OBJECT(cms. ) COUNTS PER 50 

SCAN 1 


AIR 

941 

0.0 

885 

0.5 

782 

1.0 

328 

1.5 

446 

2.0 

506 

2.5 

536 

3.0 

527 

3.5 

488 

4.0 

476 

4.5 

478 

5.0 

566 

5.5 

534 

6.0 

587 

6.5 

446 

7.0 

433 

7.5 

604 

8.0 

499 

8.5 

558 

9.0 

455 

9.5 

362 

O 

o 

765 


SECS 
SCAN 2. 
945 
846 
745 
335 
403 

525 
586 
581 
466 
451 
639 
549 
515 
554 
502 
543 

526 
520 
541 
459 
393 
812. 
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APPENDIX B 



Oject 1. 
Hollow cylinder 
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Object 2. 
Water in a beaker 


58 



59 



Object 3. 

Small cylinder in 

a large cylinder 
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n n n n n n n n n 


c 


APPENDIX C 


C PROGRAM FOR THE RECONSTRUCTION OF THE IMAGE 

C N--NUnBER' OF PROJECTIONS IN SCAN (PARALLEL) 

C M--NUMBER OF RAYS PER VI EU (FAN BEAM) 

C LITN--NUNBER OF PROJECTIONS IN SCAN(FAN BEAM) 

C LITM--NUMBER OF RAYS PER PROJECTION(PARALLEL) 

C LITH--FAN BEAM DATA (ARRAY WITH M+1 BY LITN+1 ELEMENTS) 

C D-“DISTANCE FROM SOURCE TO CENTRE OF OBJECT 

C T~~DIAMETER OF OBJECT 

C LAMMAX— MAXIMUM POSSIBLE VALUE OF LAMBDA 

C LITA--DISTANCE BETWEEN PARALLEL RAYS IN EACH VIEW 

C ALPHA---ANGULAR SPACING BETWEEN PARALLEL VI EU( RADI ANS ) 

C LITP-->PARALLEL BEAM DATA(ARRAY WITH LITM+1 BY N+1 ELEMENTS) 

C LITU-~FILTER FUNCTION DATA(ARRAy WITH LITM+1 ELEMENTS) 

C LITC--CONVOLUTED DATA(ARRAY WITH LITM+1 BY N+1 ELEMENTS) 

C LITF--BACK-PROJECTED DATA (ARRAY WITH LITH+1 ELEMENTS) 

C DELLAM--SPACING BETWEEN FAN BEAMSCIN TERMS OF LAMBDA) 

C 

REAL LITPC50,91) , LITU( 50 ) , SUM , LITC(50 , 9 1 ) , LITF ( 50 , 50 ) , T , THETA 
REAL YY,Z,Y(50),X(50),SPRIME,LITR 
INTEGER LITM,N, J ,K,L,TITLE(20) ,TTL(20) 

INTEGER F I LCBP , DATA , CONVLV .PRINT, ANS , PRT , XMAX , YMAX 
INTEGER M,LINE(50) 


FILE NUMBER VARIABLES 

0 __ 

FILCBP=1 

DATA=2 

CONVLV=3 

PRINT«4 

PIE=4 .*ATAN(1 . ) 

OPEN ( UNIT=DATA , F I LE= ' CBPPM . DAT ' ) 
READ( DATA , * )NRAY , NANGLE 


READ PARALLEL DATA 


11 F0RMAT(64A1) 

DO 40 K=l, NANGLE 

OPEN(UNIT=DATA,FILE=’PROJ21.DAT’ ) 
READ(DATA,11) (TITLE( I ) , 1=1 , 20 ) 
READ(DATA,*) (L ITP ( J , K) , J=1 ,NRAY) 
CLOSE(UNIT=DATA) 

40 CONTINUE 

C 

E=1.0 

C0NS=E 

PIXSP= 2.0*E/CNRAY-1) 
DELTH=PIE/NANGLE 
XMAX=NRAy-l 
YMAX=NRAY-1 
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(- KilHU tUUKiJbK MLlliK 

C 

OPENCUNIT=FILCBP,FILE=' FI LCBP.DAT’ ) 

READ(FILCBP,11) (TTL( I ) , I = 1 , 20 ) 

READCFILCBP,*) (LITUCK) , K=1 ,NRAY) 

DO 30 K=1,NRAY 
LITU(K)=LITU(K)/(PIXSP*PIE) 

30 CONTINUE 

CLOSE(UNIT=FILCBP) 

C . 

c 

C COMPUTE CONVOLUTION 

■C30 

DO 70 K=1,NANGLE 
DO 60 L=l,NRAy 
LITC(L,K)=0 . 0 
DO 50 J=1,NRAY 

LITC(L,K)=LITC(L,K)+(LITP(J,K)*LITUCABS(L-J)+1)) 

50 CONTINUE 

60 CONTINUE 

URITE(5,*)K 
70 CONTINUE 

c : 

C THIS PROGRAM BACK PROJECTS THE TOMOGRAPHY DATA 

C 

C CALCULATE SUPERIMPOSITION OF DATA 

C 

X(l)=-CONS 

yci)=-coNs 
DO 110 I=2,NRAY 
X(I)=X(I-1)+PIXSP 
Y(I)=YCI-1)+PIXSP 
110 CONTINUE 

SUM1=0.0 
NUMBER=0 
DO 312 K=1,NRAY 
DO 312 L=1,NRAY 
LITFCK,L)=0.0 
312 CONTINUE 

DO 120 K=1,NRAY 
DO 120 L=1,NRAY 
SUM=0.0 

DO 130 J=1,NANGLE 
AX=X(K) 

AY=YCL) 

THETA=FLOAT(J-l )*P I E/FLOAT (NANGLE) 

OTY=(AX*AX+AY*AY) 

LITR=SORT(QTY) 

IF(LITR.GT.l .0) GO TO 120 

SPRIME=(XCK)*COS(THETA)+Y(L)*SIN(THETA))/PIXSP 

G=SPRIME+FL0AT(NRAY/2+l) 

M=IFIXCFLOOR(G)) 

IF(M.GT.NRAY) LITCCM+1 , J )=0 . 0 

SUM=SUM+ ( FLOAT (M+1 ) -G) * LI TC CM, J)+(G- FLOAT (M))*LITCCM+1 , J) 

130 CONTINUE 

SUM=SUM*DELTH 

LITF(K,L)=SUM 

SUM1=SUM1+LITFCK,L) 
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120 CONTINUE 

AVE=SUM1 /NUMBER 
URITEC*,*) AVE 

C 

C THIS PROGRAM PRINTS THE TOMOGRAPHY DATA 

C 

C PRINT DATA 

C40 

YY = 99999 .0 
2=-99999.0 

DO 230 L=1,XMAX+1 
DO 220 J=1,YMAX+1 
YY=MIN(YY,LITF(L, J)) 

Z=MAX(Z,LITFCL,J)) 

220 CONTINUE 

230 CONTINUE 

C URITEC*, 35) 

C35 FORMAT ( 'NEGATIVE OF IMAGE? (0=NO , 1 = YES) ’ ) 

C ACCEPT *,ANS 

C URITEC*, 36) 

C36 FORMATC 'OPRINT IMAGE ON TERMINAL? C 0=NO , 1«YES )’ ) 

C ACCEPT *,PRT 

ANS=1 
PRT=1 

C OPEN PRINT FILE 

OPENC UNIT=PRINT, F I LE=' RESULTS. DAT’ ) 

IFCANS.EQ.l) URITECPRINT,37) 

37 FORMATC 'NEGATIVE IMAGE’) 

IFCANS.NE.l) URITECPRINT,39) 

39 FORMATC ’POSITIVE IMAGE’) 

URITECPRINT,38) CTITLECD, 1 = 1,20) 

URITECPRINT,38) (TTLCI) ,1=1,20) 

38 FORMATC 20A1) 

DO 250 L=1,XMAX+1 

DO 240 J=1,YMAX+1 

LINECO)=IFIXC0.5+99.0*CLITFCL, J)-YY)/CZ-YY)) 
IFCANS.NE.l) LINECJ)=99-LINECJ) 

240 CONTINUE 

URITECPRINT, 45) CLINECM) ,M=1 ,YMAX+1) 

45 FORMATC25I3) 

IFCPRT.EQ.l) URITEC*, 45) CLINECM) ,M=1 , YMAX+1 ) 

250 CONTINUE 

C URITECPRINT, 260) YY.Z 

260 FORMATC/ 'MINIMUM LITF<= ’, E15 . 8 , 'MAXIMUM LITF= ’ , El 5 . 8 ) 

IFCPRT.EQ.l) URITEC*, 260) YY,Z 
IFCPRT.EQ.l) URITECPRINT, 260) YY,Z 
C PRINT CENTER LINES OF PROJECTION 

URITECPRINT, 102) 

102 FORMATC 'HORIZONTAL CENTERLINE’/) 

JC=XMAX/2+l 

DO 112 I=1,XMAX+1 

URITECPRINT, 122) I , JC , LITF C I , JC ) 

122 F0RMATC1X,I3,2X,I3,5X,E15.8) 

112 CONTINUE 

URITEC*,*) CCLITFCI,J),J=JC,JC),I=1,NRAY) 

URITECPRINT, 77) AVE 
77 FORMATC/’ AVERAGE =’,E15.8) 

C50 CLOSE LISTING FILE 
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c 


c 


CLOSE (UNIT=PRINT) 

STOP 

END 


FUNCTION FLOOR(A) 


FLOOR=IFIX(A) 

IFCFLOOR.GT.A) FLOOR=FLOOR“-i . 0 

RETURN 

END 
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